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I. INTRODUCTION 



Studies of topological aspects have been attracting a great deal of attention in lattice 
QCD. Issues here include whether the index theorem relating gauge field topology and 
fermionic zero modes is well defined. Lattice fermion actions without doublers do illustrate 
the index theorem nicely. For example, one old-fashioned way to study topology in QCD is 
to use Wilson fermion [l|, |2| , although it requires fine-tuning of a mass parameter because 
of explicit breaking of chiral symmetry. The overlap fermion , which is built on the 
Wilson fermion kernel, resolves this problem since it possesses both exact chiral symmetry 
and exact zero modes related to gauge field topology, agreeing with the index theorem 



6|. The domain wall formalism [71, |8j, closely related to the overlap fermion, is also a 
useful tool to study topological effects in QCD. On the other hand, the index is hidden 
in lattice fermions with species doublers since the index effect cancels between doubling 
pairs. Although a theoretical approach to the index of the staggered fermions js ill were 



developed in j2i ll2[, it does not give an integer value from the beginning and requires a 
renormalization depending on the full ensemble of the gauge fields. Thus it is not easy to 
define the index theorem in naive, staggered and minimally-doubled fermions and thereby 
study topological effects. 



Ref. 



13| recently presented how to identify the would-be zero modes and their chiralities 



with staggered fermions away from the continuum limit using the spectral flow of a certain 
Hermitean version of the Dirac operator. The integer index obtained here correctly illustrates 
the gauge topological charge up to a factor coming from tastes. In a sequent paper jl^ . 
the author presents a new version of the overlap fermion, the "staggered overlap fermion" 
constructed from the staggered kernel. It is a two-flavor overlap fermion in four dimensions 
and the associated index of exact zero modes illustrates the gauge field topology correctly 



again. In Ref. 



15| the possibility to construct a single-flavor version of the staggered overlap 



fermion is also discussed. The key for theoretical foundation of the index theorem and the 
sequent overlap formalism here is to introduce a taste sensitive mass which is essential 
to reveal the hidden index with the spectral flow of the Hermitean operator. The new 
mass term, which assigns positive and negative masses to tastes depending on their flavor- 
chiralities, removes the cancellation of the index between tastes. When the staggered kernel 
with this kind of mass is substituted into the overlap formalism, some tastes with positive 
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mass decouple while other tastes with negative mass yield massless modes, which correspond 
to the two-flavor staggered overlap fermion. It is a universal feature for fermions with species 
doublers that you can obtain the associated overlap fermion with a proper flavored mass 
term illustrating the correct index. 

One natural extension of this approach is to apply it to the naive lattice fermion. As 
is well-known, the naive fermion includes 16 species while staggered has 4 species in four 
dimension thus it is expected we can take a parallel process to establish the integer index 
theorem without renormalization and construct overlap fermion with the naive kernel. How- 
ever one obstacle in this case is the individual species appear most easily only in momentum 
space, which leads to difficulty of identifying them as flavors and introducing flavored mass 



terms. We have a similar problem with minimally doubled fermions llLll6l-ll9ll22U23|. where 



there are two species in momentum space. So far there has been no theoretical foundation 
of the index theorem for naive and minimally doubled fermions 24 1. 

In this paper, we successfully identify species in naive and minimally doubled fermions as 
flavors by using the point splitting method, which is proposed in Ref. 20|] to deflne up and 
down quark flelds in minimally doubled fermions. By this identiflcation, we deflne proper 
flavored mass terms to extract the index in the spectral flow of the associated Hermitean 
version of the Dirac operator. Here we follow a parallel approach for identifying the would- 



13|. Then we 



be zero modes and their chiralities to that proposed in the staggered case 
flnd the spectral flow correctly illustrates the index determined by the gauge fleld topology 
in naive and minimally doubled fermions. In addition we present new versions of overlap 
fermions built from the naive fermion kernel with added flavored mass terms. Especially 
we construct a single-flavor naive overlap fermion by choosing a certain flavored mass which 
assigns negative mass to one and positive mass to all other species. One good property of this 
fermion over the minimally doubled case is it maintains hypercubic symmetry. Furthermore 
it appears to have more stability of the index against disorder than the Wilson overlap 
fermion. Thus this overlap fermion might be more suitable for practical simulations. 



II. POINT SPLITTING AND FLAVORED MASS TERMS 



With the (i-dimensional naive lattice fermion there are 2'^ species and the associated 
massless Dirac operator is anti-Hermitean. These properties are common with the staggered 
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fermion except that there are only 2*^/^ tastes in the latter case. Thus it seems easy to 
apply the approach for the index of staggered fermion to the naive fermion. In this 
spectral- flow approach the flavored mass term is essential, which gives + and — masses to 
flavors(species) depending on their chiral charges. However the species in the naive fermion 
appear naturally only in momentum space. Therefore in this case we cannot implement 
flavored mass terms straightforwardly since we need to identify species in momentum space 
as independent flavors. To resolve this problem we apply a method called "point splitting," 
originally proposed to define up and down quark fields in minimally doubled fermion 20 1. 
By using this method we obtain flavored mass terms to give proper Hermitean operators 
for the spectral flow both in minimally doubled and naive fermions. In this section we first 
introduce the point splitting for minimally doubled fermions and implement flavored mass 
terms. Secondly we perform the parallel approach for naive fermions. 



A. Minimally doubled action 



Here we introduce the point splitting method to obtain flavored mass terms in minimally 
doubled fermions. We focus only on the Karsten-Wilczek type 16|, ll7|] all through this 
paper. Now we begin with the introduction of minimally doubled fermions. The d = A 
Karsten-Wilczek action is obtained by introducing a Wilson-like term proportional to 1^4. 
Its position-space expression is 
3 



5-, 



md 



X L M=i 



+ 



sm a 



^X^^x7m {Ux 
1=1 

^(cOSa + 3)^^74^x - ^ ^^x74 {Ux,x+e^i^x+e^ + Ux,x-e^^x-eJ^ 



where the link variables satisfy Uxy = U^x- For the free theory, the associated Dirac operator 
in momentum space is given by 

3 . / 4 



Dmd{p) = i^lf^ sinp^ + cos a + 3 - ^ cospf. 



(2) 



where the parameter a adjusts the relative positions of zeros. It has only two zeros located 
a.t p = (0,0,0, ±a). These two species are not equivalent since the gamma matrices are 
differently defined between them as 7^ = r~^7^r. In the above case the transformation 
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matrix is given by F = ^7475. This means the chiral symmetry possessed by this action is 
identified flavored one given by 75 ® ra. 

The point sphtting identifles these inequivalent species as independent flavors. In this 
method each flavor fleld is deflned so that the associated propagator includes only a single 
pole. Thus we point-split the original fermion fleld by multiplying factors removing another 
undesired pole in momentum space, 

064) = ^ ( 1 + ) (3) 



2 \ sin a 

d{p + ae,) = l-T(l-'^]ij{p). (4) 



2 \ sin a 

Here u{p — ae^) and d{p + ae^) flelds correspond to the poles at p = (0, 0, 0, a) and 
(0,0,0, —a), respectively. We remark both of them yield single fermionic modes and the 
associated point-split flelds in position space are composed of the original fleld and the two 



kinds of nearest neighbors 20|. With these flavor flelds we obtain a flavor-multiplet fleld as 
following, 

\ d{p + ae^) I 

Here 75 multiplication on the original Dirac fleld is identifled as 

75^(P) j ^^(p) = (75®r3)^'(p). (6) 

Here we introduce a multiplet representation as a direct product of the Pauli matrix to 
express the 2-flavor structure. It means the flavored chiral symmetry generated by 75 (g> T3 
is exactly preserved in terms of the flavor multiplet while the flavor singlet chiral symmetry 
given by 75 ® 1 is broken in higher order terms. It is also the case of the staggered fermion 



23|. 



Now we can introduce a flavor breaking mass into minimal 



doubled fermions by a 

parallel approach to the case of staggered fermions in Refs. [13|, [ij] . We have to choose a 
flavored mass term so that it reveals the hidden index when looking into the spectral flow 
of a certain Hermitean version of the Dirac operator H{m) = 75 — M), which we will 
explain in detail in Sec. IIIII (This 75 is regarded as 75 ® T3 in terms of the flavor multiplet 
as shown in (Q.) We need to impose that 75M should be flavor-singlet as 75M = 75 Cg> 1. 
Thus the corresponding flavored mass term here is M = m (1 ® ra) with a mass parameter 
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m. To implement this flavored mass term into the action, we rewrite the mass term with 
the original Dirac fleld as 

— — siri P4. ~ 

^(p) (1 ® Ta) = uu{p — aci) — dd{p + ae^) = — ip{p)ip{p). (7) 

sm a 

It is straightforward to obtain the mass term in the position space with the link variables 
present, 

Tn — 

^md = 777—: '4'x {Ux,x+e4lpx+e4, ~ Ux.x-e^'^x-e^) ■ (8) 

2i sm a 



1 / sin^p4 



Here we also write a flavor-singlet mass term 

^{p) (1 (g) 1) ^(p) = uu{p - ae^) + dd{p + ae^) = ^ ( 1 + '^'"/^ ) ^{p)^lj{p), (9) 

I \ sm a / 

which is not required for the study of the index. The associated massive Dirac operator 
-Dmd "~ is non-Hermitean, and when gauge flelds are present the mass term does not 
commute with the kinetic term [-Dmd; ^md] 7^ 0. Thus the Dirac operator eigenvalues are 
complex. Indeed it is essential for the purpose to detect the index from the spectral flow of 
the Hermitean operator since it relies on real eigenvalues of the Dirac operator. 

In Fig. [T] we show a numerical result of complex eigenvalues of the Dirac operator for 
the d = 2 free case with a parameter a = 7r/2. Here the low-lying spectrum is split into 
two branches crossing the real axis at the magnitude of the mass parameter |m|. It means 
that the flavored mass — = diag(— m, +m) assigns — m(+m) to modes depending on 
+1(— 1) chiral charges, or equivalently +1(— 1) eigenvalues for 1 (g> T3. In other words the 
flavored mass term splits the minimally doubled fermion into two single Dirac fermions with 
— m and +m as in the staggered case. In Sec. IIIII we will see this flavored mass gives the 
spectral flow of the Hermitean operator illustrating the correct index related to the gauge 
topology. 

B. Naive action 

We can apply the same approach to the naive lattice fermion [l| to obtain a proper 
flavored mass term. The action of the naive fermions in general dimensions is simply given 

by 

1 

2^Y1 ^^^/^ {Ux,x+e^i^x+e^ - Ux,x-e^i^x-e^) ■ (10) 

x /i=l 
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FIG. 1: Complex spectrum of the free non-Hermitean Dirac operator I?md ~ -^md for the d = 2 
free field case on a 36 x 36 lattice with a mass parameter m = 1. The spectrum is split into two 
branches crossing the real axis at \m\. 



For simplicity here we consider the d = 2 naive fermions. The Dirac operator has four 
zeros (0,0), (7r,0), (0,7r) and (tt, tt) in momentum space, thus we introduce four associated 
point-split fields 









+ cospi)(l + 


COSp2)r(i)V^(p), 




-P(2)) 




- COS Pi) (1 + 


COSp2)r(2)^/'(p), 




-P{3)) 




+ cospi)(l - 


COSp2)r(3)^(p), 


^(4)(P 


-m) 


-> 


— COS Pi) (1 — 


COSp2)r(4)^(p), 



(11) 

whose locations of zeros, chiral charges and transformation matrices F giving the corre- 
sponding set of gamma matrices, 7^*'' = F|^^7^F(j), are listed in Table HI The fiavor-multiplet 
field is given by 

( ip{i){p - P{i)) \ 

i'{2){p - P{2) ) 
i^{3){P - P{3)) 
VV^(4)(P-P(4)) / 

Here the operation of 75 on the original fermion field again means the fiavored chiral 



(12) 
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label position x charge F 

1 (0,0) + 1 

2 (vr,0) - 17175 

3 (0, tt) - Z7275 

4 (7r,7r) + 75 

TABLE I: Chiral charges and transformation matrices for each of zeros in the d = 2 naive fermions 
with 71 = cji, 72 = (72 and 75 = 0-3. 

transformation in the sense of the flavor multiplet as 

/ +75 \ 

75^ (p) — > 



^(P) = (75® (r3®r3))^(p). (13) 



-75 

-75 

\ +75/ 

Here we introduce a multiplet representation as two direct products of the Pauli matrix to 
express the 4-flavor structure. Now we adopt the criterion for the proper flavored mass term 
that 75M should be flavor-singlet as 75M = 75 (g> 1. Therefore the desirable flavored mass 
for the naive fermion is given by 

#(p) (1 ® (r3 ® T^)) ^{p) = cospi cos p2ip{p)ip{p). (14) 

Please note there are also other kinds of flavored mass terms in this case: 

f(p)(l® (r3®l))*(p) = ^cospi(l + cos^p2)^(p)^(p), (15) 

f (p) (1 ® (1 ® T3)) *(p) = ^(1 + cos^pi) cosp2ipip)ip{p), (16) 

^(p) (1 ® (1 ® 1)) ^(p) = ^(1 + cos2pi)(l + cos'^p2)ip{p)ip{p). (17) 

These varieties are not useful to detect the index via the spectral flow of the Hermitean 
operator. However we will show they are essential when we search for a single-flavor naive 
overlap fermion in Sec. IIVI All of these mass terms are spread over several nearby sites, and 
therefore their position-space expressions include hopping terms with associated gauge field 
matrices. It is convenient to introduce the usual translation operators T^^ip^ = Ux,x±e^'^x±e^ 
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and = (T+^ + T_^)/2. Then the flavored mass for the index (fl^ is written in the position 
space as 

Mr,m, =mJ2 C1C2 = M„, (18) 

sym. 

where Ylsym stands for symmetric summation over the order of the factors Ci and C2, and 
m stands for a mass parameter. We can also write the extra mass terms in the position 
space as 



= ^ E (1 + ^1) (19) 

sym. 

Mi^.3 = ^E^i(l + ^2)' (20) 

sym. 

4 

sym. 

Here we comment on the possibility that the two-step hoppings in the same directions would 



^1^1 = (1 + ^1) (1+^2)- (21) 



be reduced to zero-step terms with reduction in the effects of gauge field fluctuations. 

Now let us look into the eigenvalues of the d = 2 naive Dirac operator — with 
the flavored mass in Eq. flTH]) . In Fig. H] we show a numerical result of the complex 
eigenvalues. Here the low-lying spectrum is again split into two branches crossing the real 
axis at the magnitude of the mass parameter \m\. However in this case both of the two 
branches are doubled. (We will be convinced of this doubling in Fig. |6]where it is lifted by the 
other mass terms.) It means that the flavored mass term — Mn = diag(— m, -|-m, -|-m, — m) 
assigns — m(+m) to modes depending on -|-1(— 1) chiral charges, or -|-1(— 1) eigenvalues 
for 1 (g) (ra ® T3). Thus the d = 2 naive fermion with 4 species is split into two pairs of 
Dirac fermions with — m and +m. In the next section we will show the spectral flow of the 
associated Hermitean operator gives the correct index related to the gauge topology. 

In the d = 4 case we can apply the same approach to obtain a proper flavored mass term. 
In this case there are more possibilities for flavored mass. So we will not discuss details 
here, and just show the desirable flavored mass term to reveal the index. It is given in the 
momentum space by 

))*= (^Hcosp^j^^, (22) 



and in the position space by 



Mp = mp^J]C^. (23) 

sym. 11=1 
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FIG. 2: Complex spectra of non-Hermitean Dirac operator D^ — Mrg^rs for the d = 2 free field case 
on a 36 X 36 lattice with mass parameter m = l.The spectrum is split into two doubled branches 
crossing the real axis at \m\. 

Appendix [A] is devoted to details of the flavored mass terms in c? = 4 naive fermions. In 
Fig. [H](a) for the eigenvalues of the Dirac operator we can verify the modes are split into two 
8-degenerate branches, which means the flavored mass decomposes the 4d naive fermions 
into two pairs with negative and positive masses. 



In this section we obtain the integer index related with gauge field topology for the naive 
and minimally doubled fermions. As in the cases of Wilson Isl and staggered fermions 

nn 

[131 . Il4j | we utilize the spectral flow of the Hermitean operators with the flavored mass terms 
introduced in Sec. [Tll Here the would-be zero modes of the Dirac operators are identified as 
low-lying crossings of the eigenvalue flow of the Hermitean operators. 

Let us begin with explaining what the spectral flow is. In the continuum field theory 
the index is defined as the difference between the numbers of zero modes of the massless 
Dirac operator with positive and negative chirality, ?7,+ and ?t,_. The statement of the index 
theorem is that the index is just equal to a topological charge Q of a background gauge 
configuration up to a sign factor depending on its dimensionality, 



Here the question is how to obtain the index of the Dirac operator. We can of course 



III. SPECTRAL FLOW AND THE INDEX THEOREM 




(24) 
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calculate the zero-mode chiralities straightforwardly, but there is a useful way called spectral 
flow. To introduce it we first introduce a certain Hermitean version of the Dirac operator 

H{m) =-f5{D -m). (25) 

where any zero modes of the Dirac operator with ± chirality correspond to some eigenmodes 
of this Hermitean operator with eigenvalues A(m) = =pm. If we now consider the flow of the 
eigenvalues A(m) as the mass varies, those corresponding to zero modes will cross the origin 
with slopes =f1 depending on their ± chirality. The non-zero eigenmodes of D, in contrast, 
occur in pairs which are mixed by H and cannot cross zero. Therefore the index of the Dirac 
operator is given by minus the spectral flow of the Hermitean operator, which stands for the 
net number of eigenvalues crossing the origin, counted with sign ± depending on the slope. 

The index with lattice Wilson fermions p] can similarly be obtained from the spectral 
flow, which in this case means the net number of eigenvalues crossing zero at low-lying values 
of m, counted with signs of the slopes. In the continuum limit, we are only interested in the 
crossings at small mass; the massive doublers also eventually cross zero, but only for large 
values of m. However for lattice fermions with massless species doublers the index cancels 
between pairs, so that an eigenvalue flow with a simple mass term does not properly capture 

resolved by introducing a proper 



3,Q 



The associated Hermitean 



gauge field topology. This difficulty has been recent 
flavored ("tasted") mass term in staggered fermions 
operator is given by 

H^tim) = Tss (At - "^rssFs) = T^^D^t - mT^, (26) 

where Fss and Fs stand for Fss = 75 ® 75 and 75 ® 1 with discretization error 0{a) in 
terms of a taste multiplet. Indeed it was shown that the index of the staggered fermion is 
given by minus the spectral flow of this Hermitean operator. And it was also pointed out 
that it illustrates correctly the gauge topology up to a integer factor from the number of 
tastes as Index(Dst) = 2'^/^(— 1)"'/^(5. Here the spectral flow again means the net number 
of eigenvalues crossing zero at low-lying values of m, counted with signs of the slopes. 
Thus the theoretical foundation of the index theorem with staggered fermions is established 
without a renormalization depending on the gauge ensemble. Here the mass part in the 
Hermitean operator Eq. (126|) is approximately proportional to flavor-singlet gamma-5. Note 
we constructed flavored mass terms for our cases in the previous section following this. 
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Now we can symbolically write a formula for the index as 

Index(L)) = -Spectral fiow(if). (27) 

It is quite natural to consider whether this formula is also available to detect the index of 
minimally doubled and naive fermions with the flavored mass terms we proposed. We will 
from now show this spectral flow method can be also applied to these cases. The associated 
Hermitean operators for minimally doubled and naive fermions are given by 

Hn,d{m) = 75 (^md - M^d), (28) 

Hnim) = -f,iD^ - M„), (29) 

where the matrix 75 is regarded as a flavored one, 75®r3 for minimally doubled fermions and 
75 ® (t3 eg) T3) for two dimensional naive fermions in terms of the flavor multiplet as shown in 
Sec. ini The desirable flavored mass terms and for (i = 2 have been already given 
in Eq. ([8]) and f lTSj) . For now we focus on the two dimensional case. 

We then numerically calculate the eigenvalue flows of two dimensional minimally doubled 
and naive fermions. We consider background configurations proposed in ^ for the staggered 
case [l3|: we start with a smooth U{1) gauge field with topological charge Q, 



, 1 (X2 = 1,2,--- ,L- 1) 

TT _ ZUJX2 TT _ } \ Z , , , I , . 

where L is the lattice size and u is the curvature given hj u = 2'kQ. Then, to emulate 
a typical gauge configuration of a practical simulation, we introduce disorder effects to 
link variables by random phase factors, Ux,y e^^'^'^Ux^y, where r^^y is a random number 
uniformly distributed in [—Sir, 611] . The parameter S determines the magnitude of disorder. 

Fig. El^a) shows the eigenvalue flow of the minimally doubled Hermitean operator fl2Sl) . 
It is calculated with a Q = 1 and 6 = 0.25 background configuration on a 16 x 16 lattice. 
There are two low-lying crossings around m = with positive slopes, which correspond to 
would-be zero modes. With the formula fl27j) . it means the index of the Dirac operator of 
the minimally doubled fermion in this case is —2. This result is consistent with the index 
theorem for the minimally doubled fermions given by 

Index(D^d) = 2(-l)'^/2g, (31) 
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FIG. 3: Spectral flows of (a) Minimally doubled and (b) naive Hermitean operators with a Q = 1, 
6 = 0.25 background configuration on a 16 x 16 lattice. Two single crossings with positive slopes 
are seen in (a), which means the index is —2. Two doubled crossings with positive slopes are seen 
in (b), which means the index is —4. 



(a) (b) 




-3-2-10123 -3-2-10123 



Mass Mass 



FIG. 4: Spectral flows of (a) Minimally doubled and (b) naive Hermitean operators with a, Q = 2, 
6 = 0.2 background configuration on a 16 x 16 lattice. Six single crossings with positive slopes 
and two single crossings with negative slopes are seen in (a), which means the index is —4. Six 
doubled crossings with positive slopes and two doubled crossings with negative slopes are seen in 
(b), which means the index is —8. 

which contains a factor 2 reflecting two species. This relation is also satisfled by cases with 
other topological charges, as shown in Fig. HJ^a) for the case for Q = 2. Here the net number 
of crossings counted with ± depending on the slopes is 4. It means the corresponding index 
is —4, which is consistent with ( 13T|) . We also emphasize that there is a clear separation 
between low- and high-lying crossings in Fig. where low-lying ones are localized about 
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m = and high-lying ones are located at large |m|. Although further numerical study with 
realistic gauge configurations is required, it indicates the zero modes and the index tend 
to be robust against randomness of the gauge configuration. Now we have established the 
index theorem with minimally doubled fermions. 

Next results for the naive fermion case are shown in Fig. [3^b). The calculation is done 
with the same background configuration as the minimally doubled case. Fig. [3]^b) shows 
the eigenvalue fiow of the naive Hermitean operator fl29|) . There are two doubled crossings 
around the origin: Here we can verify they are doubled and there totally exist four crossings 
as shown in Fig. [Tl^a) by introducing other kinds of mass terms given in Sec. [Tll Again with 
the formula ( 1271) we obtain the index of the Dirac operator of the naive fermion in this case 
is —4. This result satisfies the index theorem for the naive fermion given by 

lndex{D^) = (32) 

where it contains a factor 2*^ refiecting 2*^ species. This theorem is also satisfied by the 
cases with other topological charges, as shown in Fig. IH^b) for Q = 2. Here the spectral 
fiow is 8 since all the crossings are doubled in this case. It means the corresponding index 
is —8, which is consistent with fl32|) . There is also large separation between low-lying and 
high-lying crossings. It indicates the zero modes and the index tend to be robust as in the 
case of minimally doubled fermions. Now we have established the index theorem with the 
naive fermion. 

We remark both of the Hermitean operators, (128|) and (!29|) . satisfy a relation of 
75if(m)75 = —H{—m), refiecting the fact that the original actions have an exact chiral 
symmetry while the Wilson fermion does not. We can see low- lying and high-lying crossings 
get closer as the randomness effect becomes larger in the both cases, but this original chiral 
symmetry actually improves the lower bound of H'^{m) satisfying H'^ijn) = D'^ D + m?, which 
enhances the stability of the index against disorder. Thus overlap formalisms with origi- 
nally chirally symmetric fermions would be more applicable to practical simulations as we 
will discuss later. In the end of this section, let us comment on li = 4 cases. We can perform 
the same argument for the 4-dimensional minimally doubled fermions and naive fermions. 
Especially in the case of the 4d naive fermion, we can show the integer index is obtained 
from the spectral flow and it correctly illustrates the topological charge in accordance with 
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Eq. (132|) by introducing the flavored mass term 

4 

Mp = mp^J]C^, (33) 

sym. fi=l 

In Appendix |X] we will explain this mass term in detail and show the eigenvalues of the 
associated Dirac operator in Fig. E^a). 



IV. OVERLAP FORMALISM 



In this section we discuss new versions of overlap fermions constructed from the naive 
and minimally doubled Dirac kernels with the flavored mass terms. The main result here is 
we obtain a single-flavor overlap fermion with hypercubic symmetry from the naive fermion 
kernel, which may be somewhat simpler than the Wilson overlap fermion. 

Firstly we show the index of exact zero modes of the naive and minimally doubled over- 
lap fermions also illustrate the topological charge correctly. We now introduce minimally 
doubled and naive versions of overlap Dirac operators, 

Ano=l + 75^^, A.o = l+75-^, (34) 

V^md V^n 

where these Hermitean operators have been introduced in fl28p and (1291) . These depend on 
the mass parameter, which should be chosen in an appropriate region. Then we can obtain 
the corresponding Ginsparg- Wilson relations 

{75, Alio} = -Dmo75^mo, {75, Ao} = Ao75 Ao, (35) 

where, as we have discussed, the 75 here can be identifled as a flavored one in terms of 
the flavor multiplet, for example 75 — )■ 75 ® T3 for the minimally doubled fermion and 
75 75 ® ("^3 ® T3) for the d = 2 naive fermion. Note that only half of the original 
flavors(species) with negative mass are converted into physical massless modes in these 
overlap fermions while the others with positive mass become massive and decouple in the 
continuum limit. This is because the flavored mass terms we introduced assign negative 
mass to half of species and positive mass to the others. Thus there is only a single physical 
mode for the minimally doubled overlap fermion in any dimensions while there are 2*^/2 
physical modes in (i-dimensional naive overlap fermions. Here the 2*^/2 massless modes in 
the naive overlap fermion have the same chiral charge, or equivalently the same eigenvalue 
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for the "flavor- chirality" matrix given by 1 ® (ts (8x3) in the d = 2 case. As a consequence, 
the flavored Ginsparg- Wilson relations in ( l35l) reduce to unflavored relations for massless 

n 

modes as in the staggered case |14l ]. 

This reduction of flavored degrees also affects the index of the Dirac operators. We can 
obtain the indices of the minimally doubled and naive overlap Dirac operators from the 
Ginsparg- Wilson relations (135|) . 




Index(D^o) = --Tr , Index(Dno) = --Tr -j^ . (36) 




Now we can easily calculate these quantities analytically: The main part of the above equa- 
tions is a sign function as H/ \J H'^ = sgn(iJ). Thus the trace of this operator Tr (sgn(iJ)) 
gives the difference between the number of positive and negative eigenvalues at some value 
of the mass parameter. It is essential to fix the mass parameter between low- and high-lying 
crossings in the eigenvalue flows of -ffmd and H^- Then the index of the above overlap Dirac 
operator becomes just a half of that of the original Dirac operator, 

Index (Drao) = ^Index (Dmd) , Index (D^^) = ^Index (Dn) • (37) 

This relation relies on the property of the Hermitean operator '~f^H{m)'j^ = —Hl—m), as 
discussed in the end of Sec. IIIII This is also the case with the staggered overlap case. 

As shown here, we can construct a single-flavor overlap fermion from minimally doubled 
fannon., wMch howeve. lacks enough d.c.ete symmetry fa a ..nple contmuum li„.it Q 
Il9j . On the other hand, the naive overlap fermion includes lots of flavors (8 for d = 4) 
although it possess hypercubic symmetry. Now we will proceed to obtain a single-flavor 
overlap fermion from the naive fermion kernel by choosing a proper mass term. The key is 
to utilize the extra flavored mass terms given in Eqs. fITI?]) fl2D]) . 

For the purpose of constructing a single flavor naive overlap fermion with the hypercubic 
symmetry, we make use of the fact that only species (flavors) with negative mass lead to 
massless modes in the overlap formalism. It means if we introduce a flavored mass term to 
assign negative mass to one flavor and positive or zero mass to other flavors we can derive 
a single-flavor naive overlap fermion. For example, in the case of the 2d naive fermion with 
the original flavored mass fllSp . there are two degenerate species with negative mass. Now 
we need to lift this degeneracy by adding other kinds of flavored mass terms. From now we 
focus on the 2d case for a while. Now we have two possibilities of flavored mass terms to lift 




FIG. 5: (Color online) Hopping terms included in naive flavored mass terms: (a) isotropic ones, 
Mr-30r3 (blue) and Mi,g,i (red), and (b) anisotropic ones, M^-gig,! (orange) and Miig^rs (green). 

the degeneracy 

= ^ E (1 + ^i) (38) 

syrn. 
sym. 

which have been aheady shown in Eqs. ( fT9|) ( l20l) . Here we do not consider the mass term of 
Mi^i because it just gives an overall shift of the spectrum. 

There is a criterion about how to add these mass terms to the original one M^-g^T-s dH]) 
with preserving the hypercubic symmetry. As shown in Fig. Mjo), the extra mass terms 
Mi^^gand M^g^i stand for anisotropic hoppings in position space. And the only way to 
hold hypercubic symmetry is add them with the same coefficient to the original one Mr-^^r^, 
which has isotropic hoppings as shown in Fig. [St^a). Thus the possible form of the flavored 
mass to give single negative mass to the naive fermion is given by 

- Mn(c) = -Mr,(Srs " C (A'/.g^l + Mi^.g) . (40) 

where c is an overall coefficient of the two extra mass terms and we take account of the 
convention D — M . Then we find the simplest case c = 1 leads to single negative mass as 
following, 

- Mn(c = 1) = diag(-3, 1, 1, 1) = -Mn. (41) 
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FIG. 6: Complex spectra of the naive Dirac operators for the d = 2 free field case on a 36 x 36 
lattice: (a) c = 0.2, (b) c = 0.5 and (c) c = 1. In (a) the doubled negative branch is lifted and one 
of them goes to the positive direction. In (b) this branch enters a positive range. In (c) it coincides 
with the positive branch. 

This assigns only one negative mass to the eigenvalues of the Dirac operator D^^ — M^. Figure 
M^c) shows there is a single negative branch and a tripled positive branch. This means there 
is only one flavor with negative mass. We also depict figures for c = 0.2 in Fig. M^a) and 
c = 0.5 in Fig. [6](b) to convince you that one of the branch is singled out and the others are 
tripled. It describes the situation that the doubled negative branch is split and one branch 
goes towards the the positive branch. Then for the case of c ~ 0.5 this branch enters a 
positive range. Thus it is clear that the negative branch with c = 1 is singled out and the 
positive one should be triply degenerate. 

Now we find one example of the single-flavor naive overlap operator as following, 

, (42) 



(C„- rA/„)t(Z)_, _ rAf„) 

where we introduce an overall mass parameter r. Here we also study the spectral flow of the 
Hermitean operator with this mass term. Indeed it gives us a consistency check between the 
number of flavors and the index: Now we have a single flavor overlap fermion f H2|) . thus the 
associated index should be exactly equal to minus the topological charge without a flavor 
factor in this dimension. The eigenvalue flow of the naive Hermitean operator with c = 0.2, 
c = 0.5 and c = 1 for Q = 1 are depicted in Fig. [7l It is obvious that the doubled flows 
are separated in the case of c = 0.2, and we find two of the four flows no longer cross zeros 
for c = 0.5. Thus only two crossings remain in the case of c = 1. As we have discussed in 
Eq. (!36|) . the index of the overlap version is given by minus half difference of positive and 
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FIG. 7: Spectral flows of the d = 2 naive Hermitean operators with aQ = l, 5 = 0.25 background 
configuration on a 16 x 16 lattice: (a) c = 0.2, (b) c = 0.5 and (c) c = 1. The doubled flows are 
lifted in (a). Two of the four flows no longer cross zeros in (b). The two single crossings are shown 
in (c). 

negative eigenvalues of the Hermitean version of the original Dirac operator at some mass 
parameter. It means if you take an overall mass parameter r in Eq. (H2l) between low and 
high-lying crossings, the index of exact zero modes of the d = 2 naive overlap fermion ( 142|1 
with c = 1 is given by —1 for Q = 1- It generally means 

Index (^no) = -g, (43) 

which is a quarter of the original index of the 4-species naive fermion and there is no flavor 
factor. We can also check this formula for other topological charges. Thus we checked the 
consistency with the index. In the extension to general dimensions the theorem is given by 

Index (^no) = ^Index (D,) = {-ly/^Q. (44) 

Now we can also follow same approach in 4 dimensions, which is slightly more complicated. 
Here we do not show details but just present the flavored mass to give only one negative 
mass with preserving hypercubic symmetry, 

Mp + Mt + My + Ma, (45) 



20 



where the meaning of each mass term is shown in Appendix. |Xl In Fig. M^c) the eigenval- 
ues of the Dirac operator with this mass term are depicted. Here the negative branch is 
non-degenerate and the positive one is 15-fold degenerate. The associated overlap fermion 
produces only one massless mode. Thus the associated overlap form is a 4d single-flavored 
naive overlap fermion with hypercubic symmetry. You can also verify that this 4d version 
also satisfies the theorem in Eq. (H4|) . 

In the end of this section, let us comment on the symmetries and practical issues of 
this fermion. As long as the hopping symmetry of the mass terms is concerned in Fig. |5l 
hypercubic symmetry is maintained. Thus, unlike the minimally doubled case, we do not 
need to fine-tune additional parameters for the continuum limit. Regarding the numerical 
speed, it should be computationally comparable to the Wilson overlap. Furthermore the 
naive overlap fermion we proposed may have stability of the index against disorder due to 
the original chiral symmetry, although we should also take account of large fluctuations 
and renormalization effects in the realistic gauge configurations due to non-nearest hopping 
terms 26|. Now we have three examples for overlap varieties: Wilson, Staggered and Naive. 
Which is the best among them is an open question. 



V. SUMMARY AND DISCUSSION 



In this paper we have shown how the index theorem is realized in naive and minimally 
doubled fermions by considering the spectral flow of the Hermitean version of Dirac op- 
erators. The key is to make use of a point splitting for flavored mass terms. We also 
presented a new version of overlap fermions composed from the naive fermion kernel, which 
is single-flavored and maintains the hypercubic symmetry essential for a good continuum 
limit. 

In Sec. [Til we introduce the point-splitting method to identify species in momentum space 
as flavors. By using this method, we succeed to define the proper flavored mass terms in min- 
imally doubled and naive fermions. In Sec. IIIII we study the spectral flow of the Hermitean 
operators in these fermions. Then it is shown that the spectral flow correctly illustrates 
the integer index determined by gauge field topology both in naive and minimally doubled 
fermions up to overall integer factors reflecting the number of species. In Sec. IIVI we study 
new versions of overlap fermions composed from these fermion kernels with flavored mass 



21 



terms. These fermions satisfy Ginsparg- Wilson relations instead of usual chiral symmetry. 
We show the topological charge can be also obtained as the index of exact zero modes of 
the naive and minimally doubled overlap fermions. Then we show by choosing the flavored 
mass term so that only one of species has negative mass, the associated overlap fermion 
produces one massless fermionic mode. This single-flavor naive overlap fermion should pos- 
sess the hypercubic symmetry and other properties possessed by the original naive fermion. 
Furthermore it may possess more stability for the index against disorder than Wilson over- 
lap fermion since the naive fermion originally has chiral symmetry. However we also need 
to consider potential large fluctuations and renormalization effects in the realistic gauge 
configurations. Thus further investigation is required to show whether this type of overlap 
fermions has advantages for practical simulations over the Wilson overlap fermion. 

In the end of this paper, let us discuss realistic possibility to use this "Naive Overlap 
fermion" . Compared to Wilson overlap fermion, this one begins with only a simple kinetic 
term, or just a "naive action" . Once you fix a flavored mass term such as the one we have 
shown, to simulate QCD with this fermion may be easier than with the Wilson overlap. On 
the other hand, the staggered overlap fermion is based on one-component fermionic filed 
per site while our formulation uses four-component fields. This means staggered overlap 
fermion is likely to be somewhat better than the naive overlap in terms of numerical expense. 
However, to fully answer to this question we need more detailed numerical research. 
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Appendix A: Naive flavored mass in 4 dimension 



In this appendix we show details of the point sphtting and flavored mass terms for the 
d = 4 naive fermion. We introduce 16 point-spht flelds, corresponding to 16 species doublers 
of the d = 4 naive fermions, 

^(1)(P-P(1)) = ^(1 + COSpi)(l + COSp3)(l + C0SP3)(1 + COSp4)r(i)^(p), 
^(2)(P-P(2)) = ^(1 -COSpi)(l +C0SP2)(1 +C0SP3)(1 + COSp4)r(2)V^(p), 
V^(3)(P-P(3)) = ^(l + COSpi)(l-COSp2)(l+COSp3)(l + COSp4)r(3)^/'(p), 

^(16)(P-P(16)) = ^(1 - COSpi)(l - C0SP2)(1 - COSp3)(l - COSp4)r(i6)^(p), (Al) 

where the positions of zeros in the momentum space, chiral charges and definitions of trans- 
formation matrices r(j) are hsted in Table Ull Here a set of gamma matrices 7)^*^ defined for 
each zero is given by this r(j) as r^|^7^r(j) = 7^*''. We classify these zeros depending on this 
r(i). For example we denote A: axial for the case of r(2) = ^7175 while we assign T: Tensor 
for r(4) = 27172- We introduce a fiavor multiplet field with 16 components as 



P{i)) 



\ 



(A2) 



\^(16)(P 



P(16)) / 
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label 


position 


X charge 


r 


type 


1 


(0,0,0,0) 


+ 


1 


S 


2 


(vr, 0,0,0) 




11-10 


A 


3 


(0,7r,0,0) 




i ^ 10 


A 


4 


(7r,7r,0,0) 


+ 


ZTl To 


T 


5 


(0,0,7r,0) 




" /o /o 


A 


6 


(7r,0,7r,0) 


+ 


" / i /o 


T 


7 


(0,7r,7r,0) 


+ 


iT9T'? 


T 


8 


fvr, vr, vr, 0) 






V 


9 


(0, 0,0, vr) 




ZT/iTc; 


A 


10 


(vr,0,0,7r) 


+ 


iTi T/1 
" /I /4 


T 


11 


fO, 7r,0, vr) 


+ 




T 


12 


(tt TT vr) 




Tq 
;o 


V 


13 


(0,0,vr,vr) 


+ 


^73 74 


T 


14 


(vr, 0, vr, vr) 




72 


V 


15 


(0,vr,vr,vr) 




71 


V 


16 


(vr, vr, vr, vr) 


+ 


75 


P 



TABLE II: Positions of zeros, chiral charges and definitions of transformation matrices for the 
d = A naive fermions. Letters of S, V, T, V and P stand for Scalar, Vector, Tensor, Axial-vector 
and Pseudo-scalar, respectively. 
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Then 16 flavored mass terms are given by 



S: f (1 ® (1 1 (g) 1 ® 1)) * = ^ ^JJ (1 + cos2p^)j ?/;^ 

V: ^(1® (r3®l(g)l®l))^ = ^cospi ^JJ (l + cos^|>;,) j V^^/' 

T: (rg^Tg® 1® 1))* = ^cospicospa + ^°*^^^/^) j 

A : (1 ® (1 ® Ts ® Tg ® ra)) ^ = ^(1 + cos^pi) | fl cosp^ J ^j^j 



etc. 



P : ^ (1 ® (rg ® Tg ® Tg ® Tg)) ^ = ^ JJ COSp^ ^ ?/^?/^ 

where we introduce a multiplet representation as four direct products of the Pauh matrix to 
express the 16 flavor structure. In this representation the chiral transformation matrix 75 
is converted to 75 ® (rg (>>> T3 ® T3 ® T3), under which the naive action is invariant. Here the 
V-type includes 4 varieties while the T-type has 6 varieties and the A-type 4 varieties, giving 
then the total number as 16. Each of these varieties breaks the hypercubic symmetry of the 
lattice since their hoppings are anisotropic. Thus we need to take a proper combination of 
them in order to construct a flavored mass term with hypercubic symmetry. For example 
we consider the sum of the 4 varieties with the same ratios for the V-type as in the case of 
d = 2. It is also the case with the T- and A-types. Thus the 5 types of flavored masses with 
the symmetry in terms of the original fermion field are given by 

^s=5En(i+^3' (A3) 

sym. n=l 
perm. sym. v(^^i) 

= ^ E E ^M^^ (1 + ^1) (1 + ' 

perm. sym. 

perm. sym. 

4 

Mp = mp^J]C^, (AT) 

sym. /i=l 
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(a) (b) (c) 




Real part Real part Real part 



FIG. 8: Complex spectra of non-Hermitean Dirac operators for the d = A free field case in mo- 
mentum space with 16^ grids of the brillouin zone, (a) — Mp. (b) — (Mp + O.IMa). (c) 
L>n - (Mp + Mv + Mt + Ma). 

where J2perm Kieans summation over permutations of the space-time indices. 

Now we derive the flavored mass terms required to detect the index from the spectral 
flow of the Hermitean operator. As in the d = 2 case, it should be constructed so that the 
associated Hermitean operator has a flavor-singlet mass part as 75M ~ 75 (g> (1 CS> 1 1 1). 
Such a mass term is just the P-type mass ( 1A7I) . Thus the flavored mass term for the 
Hermitean operator is given by 

4 

Mp = mp5^J]C^. (A8) 

sym. /^=1 

With the Hermitean operator = 75 (-Dn — Mp), we reveal the index theorem with the 
naive fermion as in the d = 2 case. Here we only show the figure for eigenvalues of the free 
Dirac operator — Mp in Fig. IH!^a). The mass term splits the modes into two branches, 
which are 8 fold degenerate. If we introduce other types of mass terms, the degeneracy is 
lifted as seen in Fig. Mjo)- 

Next we show the flavored mass term to yield a single-flavor naive overlap fermion in 4d. 
As in the case of 2d there are some possibilities to realize it. The simplest example of the 
mass term to yield a single-flavor naive overlap fermion with hypercubic symmetry is given 
by 

Mp + Mv + Mt + Ma. (A9) 
The eigenvalues of the Dirac operator with this mass term is depicted in Fig. [8](c). Here 
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the negative branch is single and the positive one is 15 fold degenerate. It means the 
associated overlap fermion produces only one massless mode. Furthermore, as seen from 
the construction, it should possess hypercubic symmetry. Thus it could be a useful fermion 
action for lattice QCD simulation. 

We now comment on some connections with other works. The P-type mass term actually 
corresponds to the modified gamma matrix = (nu=i cosp^)75, which is introduced in 



Ref. 



25| to study the effects of topological charge with the 4d naive fermion. In some sense 



we have given its theoretical foundation by the point-splitting method. It is also similar to 
Ts ~ 75 ® 1 of the staggered fermions [23|. Regarding the T-type masses, they are similar 
to the flavored mass terms for the single-flavor staggered overlap fermion proposed in 15|. 
These similarities are reasonable since the staggered fermions can be obtained from the naive 
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